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一 种 非 均 匀 图 滤波 器 组 的 设计 方法 
卢 军 志 ， 薪 俊 正 


(桂林 电子 科技 大 学 信息 与 通信 学 院 , 广 西 桂林 541004) 


摘 要 :针对 图 频率 分 布 不 均匀 的 情况 ,提出 一 种 非 均 匀 图 滤波 器 组 的 设计 方法 。 根 据 图 频率 的 分 布 特性 ,设计 了 具备 良 
好 频率 选择 特性 及 项 点 域 局 部 特性 的 非 均匀 分 析 滤 波 器 。 先 用 低 阶 非 多 项 式 滤 波 融 近似 高 阶 多 项 式 滤 波 咒 ,再 在 给 定 分 


析 滤 波 器 和 子 带 信号 的 前 提 下 ,将 重 构 问题 归结 成 一 个 最 小 二 乘 问题 。 由 于 直接 求解 优化 问题 会 涉及 和 矩阵 求 逆 , 为 了 解决 
在 大 规模 图 的 情况 下 计算 复杂 度 高 的 问题 ,采用 一 种 预 处 理 梯度 法 对 优化 问题 进行 迭代 求解 , 且 能 分 布 式 实 现 。 仿 真 结 
表明 ,通过 该 方法 设计 的 非 均 匀 图 滤波 需 组 能 实现 完全 重 构 , 上 且 其 分 析 滤波 器 具备 良好 频率 选择 特性 及 顶点 域 局 部 特性 。 
与 其 他 迭代 方法 相 比 ,该 方法 具有 和 较 快 的 收敛 速度 和 较 低 的 计算 成 本 。 
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A design method of non-uniform graph filter banks 
LU Junzhi, JIANG Junzheng 


(School of Information and Communication, Guilin University of Electronic Technology, Guilin 541004, China) 


sn the non-polynomial form filter with low order by approximating the polynomial form with high order. Secondly, given 


non-uniform analysis filter and sub-band signal, the reconstruction problem could be formulated a least square problem. 
avoid the high calculation cost of matrix inverse in this optimization when graph is with large scale, a precondition gradi- 
descent method is proposed to solve this problem and which can be implemented in distributed manner. Numerical re- 
Tlts show that the non-uniform graph filter banks proposed in this paper can achieve perfect reconstruction and have well 
(séquency selection characteristic and localized property in vertex domain. Compared with existing iteration methods, the 
proposed algorithm has the faster convergence rate and lower calculation cost. 

Key words: precondition gradient method; graph filter banks; non-uniform graph filter; optimization problem; distributed it- 


eration algorithm 


近年 来 , 随 着 科技 和 社会 的 不 断 发 展 ,信息 数据 
呈现 出 海量 化 多样 化 、 非 规则 化 等 特点 。 传 统 的 信 
号 处 理 很 难 对 非 规 则 信号 进行 分 析 处 理 ,即使 可 以 ， 
往往 也 会 忽略 网 络 本 身 的 一 些 结构 特性 。 为 了 更 好 
地 处 理 非 规则 信号 ,研究 者 结合 图 论 等 知识 提出 了 图 
信号 处 理 , 它 能 结合 数据 本 身 的 结构 特性 对 数据 进行 
分 析 , 在 处 理 无 线 传 感 带 网 络 信 号 等 非 规 则 信号 都 取 
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得 了 不 错 的 效果 "”。 作 为 传统 信号 处 理 的 一 种 拓 
展 ,图 信号 人 处理 同样 也 存在 全 里 叶 变 换 、 滤 波 器 及 滤 
波 器 组 等 概念 ,其 中 图 滤波 器 组 因 其 具备 多 分 辩 分 析 
能 力 并 能 对 信号 进行 多 速率 信号 处 理 而 备 受 研究 者 
的 关注 。 

关于 图 滤波 带 组 的 研究 是 先 从 两 通道 开始 的 。 
2006 年 一 种 可 逆 图 滤波 器 组 的 设计 方法 中 被 提出 ， 
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并 用 于 对 传感器 网 络 上 信和 号 进行 多 分 辩 分 析 和 处 理 。 
Narang 等 趾 在 2012 年 研究 了 两 通道 临界 采样 正 交 
图 滤波 器 组 的 设计 ,通过 用 切 比 雪夫 逼近 的 方法 设计 
出 低 通 分 析 滤波 器 ,根据 正 交 特 性 求解 出 其 余 滤 波 
器 。 该 方法 缺点 在 于 不 能 满足 紧 支撑 特性 。 因 此 ， 
2013 年 , Narang 等 中 对 正 交 性 进行 了 松弛 , 变 成 双 
正 交 ,使 整个 滤波 器 组 系统 满足 完全 重 构 特性 和 紧 支 
撑 特 性 ,但 其 方法 是 基于 对 半 带 滤波 器 的 分 解 ,因此 
无 法 保证 良好 的 频谱 选择 性 。 随 后 ,研究 者 们 采用 不 
同 的 方法 "站 设计 滤波 器 组 ,以 期 达到 权衡 频谱 选择 
性 和 完全 重 构 特 性 的 目的 。 渐 渐 地 研究 者 的 关注 点 
从 两 通道 转移 到 M 通道 上 。2014 年 Tanakal 提出 
了 M 通道 过 采样 图 滤波 器 组 的 设计 方法 ,通过 构造 
过 拉 采 样 普 拉 斯 矩阵 实现 过 采样 ,最 后 实现 完全 重 
构 ? 此 后 ,一 些 基于 循环 图 的 过 采样 图 滤波 器 组 的 设 
计 郁 法 "中 被 提出 。2019 年 ,一 种 M 通道 临界 采样 
图 滤波 器 组 的 设计 方法 "被 提出 ,通过 用 插值 操作 
隘 代 重 构 部 分 来 实现 重 构 。 在 目前 的 设计 中 ,无 论 是 
宙 台 道 还 是 M 通道 ,分 析 滤 波 器 组 的 频带 划分 都 是 
均等 的 ,但 由 于 图 信号 的 频率 是 离散 的 ,频率 的 分 布 
可 能 会 出 现 分 布 不 均匀 且 部 分 集中 的 情况 ,这 时 非 均 
算 交 频率 划分 会 更 有 意义 。 
忆 在 图 滤波 器 组 的 研究 中 ,分 布 式 实现 "一直 都 
是 妆 个 很 重要 的 课题 ,特别 是 将 其 应 用 到 大 规模 的 网 
络 数据 处 理 时 ,分 布 式 实现 能 够 节省 计算 成 本 。 目 前 
大 部 分 的 图 滤波 器 组 都 采用 切 比 雪夫 多 项 式 逼 近 的 
儿 中 来 设计 ,这 种 方法 通过 增 大 阶 数 获 得 更 好 的 
频率 选择 特性 ,在 阶 数 小 的 情况 下 ,频率 选择 特性 较 
差 ? 但 具备 较 好 的 稀 琉 特性 ,能 进行 分 布 式 实现 。 然 
而 ,这 种 分 布 式 仅 限于 多 项 式 形式 下 的 滤波 器 ,并 不 
适用 于 非 多 项 式 形 式 的 滤波 器 ,如 节点 变 滤波 器 
组 5 等 。2019 年 ,Jiang 等 59 提 出 了 一 种 非 下 采样 
两 通道 图 滤波 器 组 的 分 布 式 实现 方法 ,滤波 器 的 形式 
不 再 局 限于 多 项 式 形式 ,也 可 以 是 非 多 项 式 形式 ,更 
具 普 适 性 。 该 方法 将 重 构 问 题 归结 为 一 个 最 小 二 乘 
问题 ,利用 二 阶 信息 及 近似 牛顿 法 求解 ,其 中 用 和 矩阵 
近似 理论 和 图 的 结构 特性 对 海 森 矩阵 的 道 进行 近似 ， 
从 而 实现 分 布 式 。 

针对 频率 的 分 布 特性 ,提出 了 一 种 非 均匀 图 滤波 
器 组 的 设计 方法 ,并 将 其 进行 了 分 布 式 实现 。 根 据 图 
频率 分 布 不 均匀 的 情况 ,通过 优化 方法 设计 出 具备 良 
好 频率 特性 的 多 项 式 形式 的 分 析 滤 波 器 ,用 非 多 项 式 
形式 的 滤波 器 对 其 进行 近似 。 在 已 知 子 带 信号 和 分 
析 滤波 器 的 前 提 下 ,将 图 滤波 器 组 的 重 构 问 题 归结 
一 个 最 小 二 乘 问题 。 为 了 避免 直接 求解 矩阵 的 道 , 采 
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滤波 需 组 的 设计 方法 


本 呈 革 月 二 
inaAiv 吕 TF 


用 预 处 理 樟 度 下 降 法 对 优化 问题 进行 迭代 求解 。 在 
迭代 过 程 中 ,未 利用 目标 函数 的 二 阶 信息 ,方法 更 
简单 。 


1 基础 知识 


图 信号 人 处理 中 ,无 向 图 G 可 由 三 维 集合 
G=(V,E,W) 
表示 , 其 中 : V 为 无 向 图 的 所 有 项 点 的 集合 
NN = 二 |V | ;EE 为 无 向 图 的 所 有 边 的 集合 ; W W— (ws) 
i，j EV 为 无 向 图 的 权重 和 矩阵。 根据 权重 和 矩阵 可 得 
无 向 图 的 度 和 矩阵 也, 其 对 角 线 上 元 素 


d; 一 -Dw jj 
进而 可 定义 拉 普 ; 拉 斯 矩阵 了 一 D 一 W 及 归 一 化 拉 普 
拉 斯 算 阵 LL, 一 I 一 D WD”。 在 无 向 图 中 , 工 , 是 
对 称 的 ,对 其 进行 特征 分 解 , 有 L, 二 UAU' ,U 为 特征 
向 量 组 成 的 矩阵 ， A = diag([Ai ,A ,AN |) 为 对 角 


和 矩阵 ， Al SA29°** AN 为 工 ， 的 寺 征 值 ， L, 的 寺 征 值 范 
围 为 [0,2]。 无 向 图 上 信号 x 定义 为 每 个 节点 上 的 


信号 值 的 集合 , x 二 [zi1 ,zs,…,zNj]」 ,如 图 1 所 示 。 
x 的 图 傅 里 叶 变 换 定义 为 * = 二 UTx，, 其 图 傅 里 叶 首 变 
换 为 x = 二 U* 。 根据 图 健 里 叶 变 换 的 定义 及 图 2 可 
知 , 图 频 域 上 的 频率 值 即 为 图 拉 普 拉 斯 矩阵 的 特征 
值 ,是 离散 的 , 且 4 二 0.4 和 4 二 0.8 处 具有 明显 的 频 
率 间隔 。 因 此 ,需要 根据 频率 的 聚集 程度 对 频谱 进行 
划分 ,进而 设计 出 非 均 匀 图 滤波 器 组 。 


局 
A 0.9 
pp P= 0.8 
. 
和 0.7 
一 |] 0.6 
~ 0.5 
0.4 
0.3 
多 二 
入 一 六 和 日 外 0 
WY 人 、 
四 
Y 。 0.1 
图 1 节点 数 为 50 的 随机 图 
图 滤波 的 定义 是 ,图 信号 x 经 过 图 滤波 器 后 得 


到 滤波 后 的 输出 ,其 表达 式 为 
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图 2 图 1 中 信号 的 图 频 域 表 示 


K 
f= Fx = >) ok)L’ Ky 


k=0 


其 中 p(k) 为 滤波 器 第 阶 的 系数 。 

根据 图 2 的 频率 分 布 特性 , 主要 考虑 的 是 三 通道 
非 欧 匀 图 滤波 器 组 。 由 于 目前 图 上 的 采样 操作 并 不 
叭 全 ,是 没有 针对 奇数 通道 的 采样 设计 ,因此 不 讨论 
图 泪 波 器 组 上 的 采样 操作 三 通道 非 均 匀 非 下 采样 
天 刘波 器 组 吉 构 框图 如 图 3 所 示 , x 、 守 分别 为 图 滤 
滤 罗 组 的 输入 和 输出 ,，{H,,H),H,} 为 分 析 滤 波 器 
zo,z1,z2) 为 滤波 器 的 子 带 系数 , zw 二 Hx， 
[0,1,2], {Go ,G1,Gs) 为 重 构 滤波 器 组 。 该 框 
架 的 输入 输出 关系 可 表示 为 

>< =Tx=(GH,GH+G,H,)x. (1) 
为 满足 其 完全 重 构 特性 ,整个 系统 的 传递 函数 了 应 
满 定 ， 


攻 
md 


,9 T= Go Ho, + 1 CI 五 ; + 
其 中 工 为 单位 矩阵 。 


| GH,=I,， (2) 


图 3 三 通道 非 下 采样 图 滤波 器 结构 组 框图 


目前 大 部 分 的 滤波 咒 组 设计 均 和 采用 切 比 雪夫 近 
似 的 方法 同时 设计 分 析 滤 波 器 和 综合 滤波 顺 。 这 些 
方法 都 是 设计 多 项 式 形 式 的 ,设计 自由 度 不 高 , 且 在 
良好 的 频率 选择 特性 与 好 的 局 部 特性 之 间 不 能 很 好 
地 权衡 。 阶 数 越 高 ,频率 选择 性 越 好 ,但 局 部 特性 越 
差 , 越 不 能 分 布 式 实现 。 而 非 多 项 式 形式 的 滤波 顺 在 
提高 设计 自由 度 的 前 提 下 尽 可 能 使 得 频率 选择 性 更 


好 ， 同时 也 拥有 良好 的 局 部 特性 ， 因此 ， 法 用 非 多 项 
式 形式 中 的 NV 图 滤波 器 作 为 分 析 滤 波 咒 ,其 设计 步 
又 为 : 

1) 通 过 约束 滤波 器 的 通 带 平坦 性 和 阻 带 衰减 得 
到 频率 选择 性 好 的 多 项 式 形 式 的 分 析 滤波 器 : ,其 
优化 问题 可 表示 为 


b d 
min| |1—b' Vol da | 1D CA)p | ?qd)， 
9 a c 


(3) 
站 
XA*]"。 式 (3) 的 前 项 是 对 滤波 器 通 带 平坦 性 的 约束 ， 
a 0 为 滤波 融 的 通 带 截止 频率 ;后 项 是 对 阻 带 衰 减 的 
约束 , c 、d 为 滤波 器 的 阻 带 截止 频率 , a 为 可 调整 参 
数 。 确 定好 滤波 器 的 通 带 、 阻 带 截止 频率 后 ,就 能 通 

过 求解 式 (3) 的 优化 问题 得 到 图 滤波 器 的 系数 。 
2) 将 设计 好 的 多 项 式 滤 波 器 当 作 理想 滤波 器 


理 , 再 用 优化 的 方法 设计 出 阶 数 较 低 且 具 有 良好 频 
率 选 择 性 的 NV 图 滤波 器 中 。 在 该 方法 中 ,通过 最 
小 化 Frobenius 范 数 
| 过 ce 天， 

得 到 最 优 的 滤波 器 系数 , Ci 为 系数 和 矩阵 。NV 滤波 
器 的 优点 是 获得 更 小 阶 数 、 更 好 频率 选择 特性 的 滤波 
需 , 同 时 所 需 设计 的 系数 变 多 ,设计 自由 度 变 高 ,设计 
难度 也 变 大 。 

图 滤波 吉 组 的 重 构 部 分 并 不 需要 像 分 析 滤 波 需 
一 样 具 有 良好 的 频谱 选择 特性 ,只 需 满 足 完全 重 构 特 
性 即 可 。 因 此 , 重 构 部 分 的 实现 也 不 一 定 需 要 通过 重 
构 滤 波 带 来 实现 , 且 目 前 并 无 太 多 关于 非 多 项 式 形式 
图 滤波 带 组 的 设计 方法 。 在 给 定 了 分 析 滤 波 絮 组 
{Ho ,Hi ,Hi) 及 子 带 信号 (zo,z1,z2) 的 前 提 下 ,将 


重 构 问题 归结 为 一 个 优化 问题 : 
min Hox —zol;+ HX zl: 
| H;,X —z; |;, (4) 
此 最 小 二 乘 问题 的 最 优 解 为 
X ”一 百 xz， (5) 
其 中 : 


H=H,o H+Hi'Hi+H,'H;; 


z= Ve 
最 优 解 的 求解 需要 进行 矩阵 求 逆 运 算 ,而 矩阵 求 逆 的 
复杂 度 是 OCN  ) 。 若 矩 阵 维 度 很 大 , 且 求 逆 和 矩阵 无 
特 丈 性质 的 情况 下 , 求 逆 的 计算 代价 很 高 , 且 可 能 存 


在 矩阵 病态 问题 ,因此 采用 一 种 预 处 理 分 布 式 梯度 下 
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降 法 来 求解 这 个 道 滤波 问题 。 
2 设计 方法 

逆 滤 波 问题 是 一 个 线性 问题 ,可 采用 梯度 法 进行 
夫 代 求解 ,并 采用 预 处 理 技术 加 速 迭代 方法 的 收敛 速 
度 。 预 处 理 技术 的 基本 思想 是 在 原本 的 线性 问题 中 
加 入 预 处 理 算 子 ,使 之 能 等 价 转换 为 新 的 线性 问题 。 
这 么 做 的 好 处 是 能 够 获得 更 快 的 收敛 速度 ,同时 在 一 
些 条 件数 很 差 的 系统 中 ,还 能 起 到 改善 条 件数 的 作 
用 。 本 人 研究 选取 的 预 处 理 算 子 P 定义 为 
>， IHG,7)|,i€EV, 


jEB(i,wH)) 


P= diag(pi), pi= 


(6) 
其 中 : BG,w (H))=={j EV,0(i,j) 之 w (HH)) 为 节 
点 亏 的 w(H) 阶 以 内 所 有 邻居 节点 的 集合 ; w(H) 为 
珑 阵 H 的 测 地 距离 , w(H) 三 0; p(i,j) 为 从 节点 i 
节点 j 所 需 的 最 短路 径 。w( 互 ) 与 滤波 器 的 阶 数 
j 闫 ,是 阶 数 的 2 倍 。 设 置 训 一 P ?HP , 类似 归 
-人 tg 拉 普 拉 斯 矩阵 的 操作 ,对 和 矩阵 了 H 进行 归 一 化 处 
理 N 本 研究 所 选 的 预 处 理 形式 是 右 预 处 理 形式 ， 


Hy =P xz 9 py:, (7) 
将 其 与 梯度 法 结合 后 ,有 
A 
= (mt1) (Gm) TF ， (7) 一 1] /2 
= —(H —P 2 
> . 4 (8) 
< TY (m++1) Py, 


和 绎 过 变换 后 ,可 得 
SF0nTD) To0) PHTGn)— 7). (9) 
-于 在 整个 方法 中 ,避免 了 复杂 矩阵 的 求 着 问 题 。 和 撼 
阵 的 测 地 距离 wa(H) 是 滤波 器 阶 数 的 2 倍 , 而 由 于 所 
设计 的 滤波 器 阶 数 小 , 故 w(H) 很 小 。 因 此 ,在 每 次 
人 欠 代 过 程 中 ,节点 只 与 其 少量 的 邻居 进行 信息 交换 ， 
整个 过 程 可 分 布 式 实现 。 该 方法 对 每 个 节点 上 上 信 
号 进行 分 布 式 重 构 的 算法 流程 如 下 : 

1) 初 始 化 庆 (0)G) = 0,; € BG,w(H))， 
Xx (0)()=0,m=0, 


计算 p= 9 [G7 
jE€EB(i,w(H)) 

while (mm 和 M) do 

3) 计 算 


TOG)=T0n)0) Tz/p; CO— 
HO, (I)/p;, 


jE€EB(i,wH)) 
4) 将 让 ?G2) 发 给 节点 i 的 w CH) 阶 以 内 所 有 
邻居 节点 ， 
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5) 从 邻居 节点 接收 由 计算 得 到 的 
xX (m1)G)。 


end while 


将 所 有 的 UM) (i) 进行 组 合 ,可 得 到 所 有 节点 
最 终 的 重 构 信号 值 。 由 算法 流程 可 看 出 ,在 整个 迭代 
过 程 中 ,计算 复杂 度 最 高 的 操作 就 是 矩阵 与 向 量 的 相 
乘 ,此 计算 复杂 度 比 矩阵 求 逆 的 计算 复杂 度 低 了 一 个 
数量 级 。 因 此 ,该 方法 的 计算 复杂 度 小 于 直接 求解 的 
复杂 度 , 且 在 图 规模 很 大 的 情况 ,这 种 优势 会 更 明显 。 


3 仿真 实验 


首先 给 出 关于 三 通道 非 均 匀 图 滤波 器 组 的 分 析 
滤波 器 组 的 设计 及 分 析 , 然 后 给 出 该 方法 与 现 有 的 分 
布 式 迭代 方法 的 对 比 。 仿 真 所 用 的 图 信号 如 图 1 所 
示 , 所 有 仿真 均 在 相同 环境 下 进行 。 

实验 1 根据 图 2 的 频率 分 布 特性 ,通过 优化 非 
均匀 分 析 滤 波 器 组 的 系数 ,设计 出 如 图 4 所 示 的 非 均 
匀 频 谱 划 分 的 分 析 滤 波 器 组 。P Ai， 分别 表示 低 
通 、 带 通 和 高 通 非 均匀 分 析 滤 波 器 的 频率 响应 ,3 个 
滤波 器 的 阶 数 都 是 16。 从 图 4 可 看 出 , 低 通 、 高 通 分 
析 滤 波 器 的 通 带 都 比较 平坦 ,频率 选择 特性 都 很 好 ， 
高 阶 多 项 式 滤波 器 虽然 频率 选择 特性 好 ,但 局 部 特性 
很 差 , 并 不 适合 应 用 于 分 布 式 方法 中 ,因此 采用 NV 
滤波 器 中 对 设计 好 的 多 项 式 滤波 器 进行 近似 。 在 实 
验 中 ,设计 了 阶 数 依次 为 1.2 和 3 的 NV 滤波 器 ,分 
别 用 Kl、K2、K3 表示 。 但 由 于 NYV 滤波 器 并 不 存在 
频率 响应 的 说 法 ,为 了 说 明 所 得 到 的 NV 滤波 器 的 滤 
波 性 能 与 图 4 的 性 能 相近 ,将 高 频 子 带 信号 进行 比较 
来 说 明 近 似 效果 ,观察 NV 滤波 器 经 过 高 频 滤 波 后 的 
得 到 的 高 频 子 带 信号 的 频谱 图 是 否 与 多 项 式 形式 滤 
波 器 一 致 ,从 而 判断 NV 形式 的 滤波 器 是 否 达到 了 较 
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幅 


图 4 设计 得 到 的 三 通道 非 均匀 分 析 滤 波 器 组 的 幅度 响应 
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好 的 近似 效果 。 图 5 为 输入 信号 经 过 不 同形 式 的 高 
通 滤波 器 后 得 到 的 子 带 信号 频 域 表 示 。 从 图 5 可 看 
出 ,NV 形式 滤波 器 滤波 得 到 的 子 带 信号 频谱 尽管 在 
少许 频率 分 量 上 会 有 幅度 上 的 小 差别 ,但 整体 上 滤波 
效果 与 多 项 式 形式 的 效果 接近 。 因 此 ,采用 NV 形式 
滤波 器 来 近似 高 阶 多 项 式 形式 滤波 器 的 结果 是 可 行 
的 ,能 取得 同样 的 效果 。 

实验 2 对 图 信号 进行 分 布 式 重 构 的 实验 仿真 ， 
总 的 迭代 次 数 为 200 ,分析 滤波 器 基于 例 1 的 多 项 式 
分 析 滤 波 器 ,是 通过 优化 设计 出 的 非 多 项 式 形式 的 分 
析 滤 波 器 。 当 滤波 器 分 别 为 K1、K2 和 K3 时 ,不 同 
迭代 方法 的 收敛 曲线 分 别 如 图 6.7 和 图 8 所 示 。 相 
对 误差 为 ‖ 庆 Gn) 一 x 目 ,/ x |。 在 仿真 中 发 现 ， 
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(c) K2 滤 波 后 的 高 频 子 带 频 域 表示 
5 相同 信号 经 
表 1 迁 代 完成 后 本 算法 在 不 同 滤波 器 阶 数 下 的 重 构 误 差 
滤波 器 ”分布 式 梯度 法 中 分 布 式 牛顿 法 D 本 方法 
K1 3.28X10 2 1; 85X10 1.02X10 于 
K2 5 140 22 321RLIO 各 6.20X 10 32 
K3 .14X10 3 4.35X10 3 3. 40X10 3 


当 步 长 取 0.9 时 ,分 布 式 梯度 法 “能 达到 最 优 的 效 
果 。 因 此 ,在 不 同 的 阶 数 下 ,该 方法 均 采 用 相同 的 步 
长 。 根 据 图 6 一 8 可 知 ,本 方法 的 迭代 收敛 次 数 比 分 
布 式 梯度 法 要 少 很 多 ,与 分 布 式 牛顿 法 "相差 不 
大 ,在 不 使 用 二 阶 信息 的 前 提 下 ,取得 了 比 一 般 分 布 
式 梯度 法 更 快 的 收敛 速度 ,表明 预 处 理 技术 能 够 有 效 
加 快 迭 代 收 敛 的 速度 ;在 阶 数 不 断 增 大 时 ,所 有 算法 
的 收敛 次 数 都 在 减少 ,原因 在 于 阶 数 增 大 ,节点 与 周 
围 邻 居 节 点 的 通信 次 数 变 多 ,获得 的 有 用 信息 增多 ， 
通过 不 断 的 信息 交换 ,就 能 更 快 地 达到 收敛 。 表 1 为 
迭代 完成 后 本 方法 在 不 同 阶 数 条 件 下 的 重 构 误 差 。 从 
表 1 可 看 出 ,本 方法 能 够 进行 分 布 式 重 构 , 且 重 构 误 差 
比分 布 式 梯度 法 ”小 ,与 分 布 式 牛 顿 法 "差距 不 大 。 
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(d) K3 滤 波 后 的 高 频 子 带 频 域 表示 


过 不 同 滤波 器 得 到 的 高 频 子 带 信号 的 频 域 表示 


4 结束 语 


针对 频率 分 布 不 均匀 的 情况 ,提出 了 一 种 非 均匀 
非 下 采样 图 滤波 噩 组 的 设计 方法 ,并 进行 了 分 布 式 实 
现 。 根 据 图 频率 的 分 布 特性 ,设计 了 具备 良好 频率 划 
分 特性 及 稀 玻 特性 的 低 阶 非 多 项 式 分 析 滤 波 器 。 在 
已 知 分 析 滤 波 器 和子 带 信和 号 的 前 提 下 ,将 图 滤波 器 组 
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图 8 采用 K3 时 ,不 同方 法 的 收敛 曲线 

的 重 构 问 题 归结 为 一 个 最 小 二 乘 问题 ,并 采用 预 处 理 
梯度 法 对 其 进行 迭代 求解 。 仿 真 结 果 表明 ,加 入 预 处 
理 算 子 后 ,能够 加 快 收敛 速度 ,上 且 比 利用 二 阶 信息 的 


迭代 方法 拥有 更 低 的 计算 成 本 。 
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